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Abstract 

In three-dimensional computational topology, the theory of normal surfaces is a tool of great 
theoretical and practical significance. Although this theory typically leads to exponential time 
algorithms, very little is known about how these algorithms perform in "typical" scenarios, or 
how far the best known theoretical bounds are from the real worst-case scenarios. Here we 
study the combinatorial and algebraic complexity of normal surfaces from both the theoretical 
and experimental viewpoints. Theoretically, we obtain new exponential lower bounds on the 
worst-case complexities in a variety of settings that are important for practical computation. 
Experimentally, we study the worst-case and average-case complexities over a comprehensive 
body of roughly three billion input triangulations. Many of our lower bounds are the first known 
exponential lower bounds in these settings, and experimental evidence suggests that many of 
our theoretical lower bounds on worst-case growth rates may indeed be asymptotically tight. 

1 Introduction 

In three-dimensional computational topology, many important problems are solved by exponential- 
time algorithms: key examples include Haken's algorithm for recognising the unknot [13], or break- 
ing down a triangulated 3-manifold into its prime decomposition [18\ I19j. This is in contrast to two 
dimensions in which many problems are solved in polynomial time, or higher dimensions in which 
important topological problems can become undecidable [11[ I24j. 

A common feature of many such three-dimensional algorithms, and the source of both their 
solvability and their exponential running times, is their use of normal surfaces. In essence, normal 
surfaces are embedded 2-dimensional surfaces that intersect the surrounding 3-dimensional triangu- 
lation in a simple fashion. Most importantly, they describe topological features using combinatorial 
data, and are thereby well-suited for algorithmic enumeration and analysis. 

Amongst the most important normal surfaces are the vertex normal surfaces. These correspond 
to the vertices of a high-dimensional polytope (called the projective solution space), and together 
they generate the space of all possible normal surfaces within the input triangulation. Many 
topological algorithms begin by enumerating all vertex normal surfaces in the input triangulation, 
and for many problems (such as unknot recognition and prime decomposition) this enumeration is 
in fact the main bottleneck for the entire algorithm. 
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One remarkable feature of many algorithms in three-dimensional computational topology is 
that, although they have extremely large theoretical worst-case complexity bounds, they appear to 
be much easier to solve in practice than these bounds suggest. For example: 

• In 1980, Thurston asked if the Weber-Seifert dodecahedral space is Haken (the precise meaning 
of this is not important here) PQ. This long-standing question became a symbolic benchmark 
for computational topology, and was only resolved by computer proof after 30 years |10j . 
At the heart of the proof was an enumeration of all vertex normal surfaces in an n = 23- 
tetrahedron triangulation: despite a prohibitive 0(16 n x poly(n))-time enumeration algorithm 
(the best available at the time) and a best known bound of 0(3.303") vertex normal surfaces 
[5j, the enumeration ran in just 5^ hours with only 1751 vertex normal surfaces in total. 

• The problem of unknot recognition is of particular interest. Modern derivatives of Haken's 
original algorithm |19j have an exponential time complexity |16j . but there is a growing dis- 
cussion as to whether a faster algorithm might exist |12[ [T5] . Certainly unknot recognition lies 
in NP [16] . and also co-NP if we assume the generalised Riemann hypothesis [21]; moreover, 
recent algorithmic developments based on linear programming now exhibit an experimen- 
tal polynomial-time behaviour [;6J. Deciding whether unknot recognition has a worst-case 
polynomial-time solution is now a major open problem in computational topology. 

This severe gap between theory and practice is still poorly understood. There appear to be two 
causes: (i) the best theoretical complexity bounds are far from tight; (ii) "pathological" inputs that 
exhibit high-complexity behaviour are rare, with "typical" inputs often far easier to work with. 

Proving such claims mathematically remains extremely elusive. Obtaining tight complexity 
bounds requires a deep interaction between topology, normal surfaces and polytope theory, and 
it is difficult to avoid making very loose estimates in at least one of these areas. Understanding 
"typical" behaviour (such as average- or generic-case complexity) is hampered by our very limited 
understanding of random 3-manifold triangulations: even the simple task of generating a random 
3-manifold triangulation with n tetrahedra has no known sub-exponential-time solution [27] . In 
this setting, experimental work plays a crucial role in understanding the realistic performance of 
algorithms, as well as the innate difficulty of the problems that they aim to solve. 

In this paper we focus our attention on the problem of enumerating all vertex normal sur- 
faces within a given n-tetrahedron input triangulation: as mentioned earlier, this is a central 
component — and often the main bottleneck — of many algorithms in computational 3-manifold 
topology. Enumeration algorithms are still evolving [BJ [8], and they are often hand-tailored to 
a particular topological problem of interest. For this reason we do not focus on the complexity of 
any specific algorithm, but instead we study two aspects of normal surface theory that affect and 
constrain all of these algorithms: 

• Combinatorial complexity: We study the total number of vertex normal surfaces within the 
input triangulation T, which we denote by o~{T). This is our main quantity of interest. It 
yields an immediate lower bound for the time complexity of any enumeration algorithm, since 
it determines the output sizeQ Moreover, o~(T) also factors into upper bounds, since modern 
enumeration algorithms are designed to run faster in situations where o~(T) is small [BjH 

• Algebraic complexity: As detailed in Section [21 each normal surface is described by a non- 
negative integer vector in IR 7n (or in some settings, R 3n ). We investigate the maximum co- 
ordinate of this vector of any vertex normal surface within the input triangulation T, which 

1 Specifically, since each vertex normal surface can be described in 0(n) space [16| . the output size is 0(a(T) x n). 
2 The tree traversal enumeration algorithm (the current state of the art) has running time 0(4™<t(T) x poly(n)). 
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we denote by k(7"). This quantity is important for the implementation of enumeration algo- 
rithms, since it affects whether we can work with fast native machine integer types or whether 
we must fall back to significantly more expensive arbitrary-precision integer arithmetic [8]. 
Moreover, k(T) features in algorithms that extend or even avoid the enumeration problem: 

— Some algorithms, such as recognising small Seifert fibred spaces [28J, require the complete 
enumeration of not just vertex normal surfaces but a much larger "lattice" of normal 
surfaces whose size is a function of k(T). 

— Some algorithms, such as determining the crosscap number of a knot [9], avoid vertex 
enumeration entirely by solving an integer program instead; here the bounds on k(T) 
feature as coefficients in the integer program, and directly affect whether the program 
can be solved using off-the-shelf integer programming software. 

In summary, by focusing our attention on the quantities o~(T) and re(T), we learn not only 
about the behaviour of current enumeration algorithms, but also about the intrinsic limits and 
behaviour of the problem that they seek to solve. 

We approach these combinatorial and algebraic quantities <j{T) and k(T) through both theory 
and experiment. Theoretically, we construct infinite "pathological" families of triangulations in 
Section[3]that establish exponential lower bounds on the worst-case scenario for both 0"(7~) and k(T). 
Experimentally, we examine both the worst case and average case behaviour of these quantities in 
Section [H using a comprehensive census of billions of input triangulations. 

Such results are highly important for practitioners in three-dimensional computational topology, 
particularly given the exponential nature of many key algorithms. Despite this, just one preliminary 
study of this type appears in the literature [3]. This scarcity of results has two causes: 

(i) The lack of large, comprehensive censuses of both "typical" and "atypical" triangulations. 

There are many censuses of 3-manifold triangulations in the literature, but these typically 
focus on well-structured triangulations with special properties (such as minimal triangulations, 
or irreducible manifolds) . Such triangulations are often easy to work with [18] , and offer little 
insight into an algorithm's worst-case (or even average-case) behaviour. 

It is only recently that large, comprehensive bodies of census data have been developed to 
study all triangulations of a given input size [U 0]. By using such censuses for our experimental 
data, we ensure that we identify pathological cases, and also gain a clear understanding of 
how common or rare they are. 

(ii) The intense computation required to study normal surfaces with such large bodies of data. 

Normal surface enumeration algorithms have enjoyed significant advances in recent years, and 
modern algorithms now run many orders of magnitude faster than their earlier counterparts 
[U [8]. The experimental work in this paper required several years of combined CPU time, 
and without recent algorithmic advances [21 01 [6] this work would not have been possible. 

The preliminary study in [3] examines only the combinatorial complexity <r(T), and works with 
a data set of roughly 150 million triangulations of closed 3-manifolds. The study in this paper is 
significantly richer, both in scope and detail: 

• We examine the algebraic complexity n(T) in addition to the combinatorial complexity cr(T); 

• We work with a comprehensive data set of almost three billion triangulations, spanning both 
closed manifolds (which are important for algorithms such as prime decomposition) and 
bounded manifolds (which are important for knot algorithms); 
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• We also examine these quantities in "optimised" settings that arise in practical computation — 
in particular, one vertex triangulations (a common optimisation used in many topological 
algorithms), and the restricted problem of enumerating only vertex normal discs (which is 
important for unknot recognition, or testing surfaces for incompressibility) . 

Our pathological families yield the first known explicit exponential lower bounds on worst-case 
complexity for the computationally important settings of bounded triangulations, closed 1-vertex 
triangulations, and normal discs. In many settings our pathological families match the experimental 
worst-case growth rates precisely, and we conjecture that the resulting exponential bounds are in 
fact exact. 

Of related note is a result of Hass et al. [T7], who establish an exponential lower bound on 
the worst-case complexity of a triangulated disc spanned by the unknot in M 3 (this has particular 
relevance for the complexity of unknot recognition). Their result operates under stricter geometric 
constraints, and it is not yet known how it translates to the more flexible setting of normal surfaces. 

We emphasise again that our experimental data sets use exhaustive censuses of all possible 
input triangulations below a given size. This paper introduces the first such census of bounded 
3-manifold triangulations in the literature, totalling over 20 billion triangulations of size n < 9. 

We use exhaustive censuses because there is no known efficient algorithm for randomly sam- 
pling large triangulations \13\ [27] , and although there are other methods for generating random 
3-manifolds |13|. 123]. nothing is known about the bias of the resulting sample of triangulations. As 
a result, although our census is large, the triangulations it contains are all relatively small. Never- 
theless, there are strong reasons to believe that our experimental results are indicative of behaviour 
for larger inputs; we discuss this further in Section 

2 Preliminaries 

By a triangulation T, we mean a collection of n abstract tetrahedra Aj = i(0123), 1 < i < n, 
some or all of whose faces are affinely identified or "glued together" in pairs; here (0123) refers 
to the four vertices of tetrahedron Aj. As a consequence of these face gluings, many tetrahedron 
edges may become identified together; we refer to the result as a single edge of the triangulation, 
and likewise with vertices. The gluings must be in a way that no edge is identified with itself 
in reverse as a result. Moreover, each tetrahedron face must be identified with one and only one 
partner (we call these internal faces), or with nothing at all (we call these boundary faces). The 
set of boundary faces is called the boundary of the triangulation and denoted by dT. If dT = 
then T is called a closed triangulation, otherwise it is said to be bounded. Not all triangulations 
(closed or bounded) represent 3-manifolds. However, unless otherwise specified, this is true for all 
triangulations presented in this paper. 

Throughout this article, the gluings of the triangles are given by a bijection of their vertices 
i(abc) i — y j(def) where the symbol i(abc) denotes triangle (abc) from tetrahedron i and the order of 
the vertices determines the gluing. A triangulation as defined above is sometimes referred to as gen- 
eralised triangulation; these are more general and flexible than simplicial complexes. An important 
case is a 1-vertex triangulation, in which all tetrahedron vertices become identified together. 

The face pairing graph of a triangulation T is the multigraph whose nodes represent tetrahedra, 
and whose arcs represent pairs of tetrahedron faces that are glued together. A face pairing graph 
may contain loops (if two faces of the same tetrahedron are glued together), and/or multiple edges 
(if two tetrahedra are joined together along more than one face). See Figure [3] for examples. 

A properly embedded surface in T is a surface s C T with no self- intersections, and whose 
boundary lies entirely within dT. A normal surface in T is a properly embedded surface that 
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Figure 1: Normal triangles and quadrilaterals within a tetrahedron. 



meets each tetrahedron A of T in a disjoint collection of triangles and quadrilaterals, each running 
between distinct edges of A, as illustrated in Figure [TJ There are four triangle types and three 
quadrilateral types in A according to which edges they meet. Within each tetrahedron there may 
be several triangles or quadrilaterals of any given type; collectively these are referred to as normal 
pieces. The intersection of a normal piece of a tetrahedron with one of its faces is called normal 
arc; each face has three arc types according to which two edges of the face an arc meets. 

Counting the number of pieces of each type for a normal surface s gives rise to a 7-tuple per 
tetrahedron of T and hence a 7n-tuple of non- negative integers describing s as a point in M>q, called 
its normal coordinates. Such a point must satisfy a set of linear homogeneous matching equations 
(one for each arc type of each internal face). These equations are necessary but not sufficient: 
the normal coordinates must also satisfy a set of combinatorial constraints called the quadrilateral 
constraints, which we discuss further in the appendix. 

The solution set to the matching equations in R>q is a polyhedral cone (the cross-section 
polytope of this cone is also known as the projective solution space). A vertex normal surface is one 
whose normal coordinates lie on an extremal ray of this polyhedral cone and, in addition, its normal 
coordinates are minimal for all integer points on this ray. Thus, there are only finitely many such 
vertex normal surfaces; every normal surface can then be expressed as a positive rational linear 
combination of these surfaces just like every point in a polyhedral cone is a positive rational linear 
combination of points on its extremal rays. This is why, when enumerating normal surfaces in a 
triangulation, we typically just consider the finite set of vertex normal surfaces. 

3 Theoretical lower bounds 

Here we establish lower bounds for the worst-case values of o~(T) and n(T), i.e., the maximum 
possible cr(T) or k(T) for an n-tetrahedron triangulation T ■ Recall that cr(T) measures the combi- 
natorial complexity, i.e., the number of vertex normal surfaces within T, and that k(T) measures 
the algebraic complexity, i.e., the maximum coordinate of any vertex normal surface of T '■ 

Few such results are known: there are no explicit lower bounds on the worst-case k(T) in the 
literature, and the only explicit lower bound on the worst-case cr(T) is given by a family of closed 
triangulations with <r(T) G G(17 n/4 ) ~ 9(2. 03 n ) [3]. In this section, we give new exponential lower 
bounds for cr(T) and k(T) in a variety of settings that hold particular relevance for key algorithms 
in computational geometry and topology. We sketch the main constructions and results here; see 
the appendix for detailed proofs. 

3.1 Closed triangulations with many normal surfaces 

Important 3-manifold algorithms such as prime decomposition often begin by converting the input 
triangulation to a 1-vertex triangulation, whereupon the subsequent processing becomes signifi- 
cantly easier [18]. The 6(2. 03 n ) family of [3] is not of this type (each triangulation has n + 1 
vertices), which raises the question of how such bounds behave in a 1-vertex setting: 
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Theorem 1. There is a family A n , n > 1, of closed 1-vertex triangulations with n tetrahedra and 
o~(A n ) = 2 n vertex normal surfaces. 

We call these triangulations binomial triangulations, because more precisely they have (?) vertex 
normal surfaces of genus k for each k = 0, . . . ,n (whereby a(A n ) = (1) = 2 n ). We construct 
each A n from n tetrahedra Ax, ■ ■ ■ , A n in the following manner. 

We begin by folding together two faces of the same tetrahedron Aj = i(0123) by the gluing 
z(012) £(013) for i = 1, ...,n. Then we identify tetrahedra A, and A i+1 by £(123) ^ (i + 1)(230) 
(where A n+ i = Ai). These gluings identify all the vertices to a single vertex, therefore An is a 
1-vertex triangulation. See Figure [3] for a picture of the face pairing graph of A n for n = 6. It can 
be shown that each A n is a closed 1-vertex triangulation of the 3-sphere. 

To see why a{A n ) = 2 n , we observe that in each tetrahedron there are two normal "subsurfaces" 
which are compatible with any other normal surface of the triangulation. One of these adds genus to 
the overall surface and the other does not. We show that the vertex normal surfaces are precisely 
combinations of these subsurfaces, whereby the binomial coefficients and 2 n growth rate easily 
follow. 

Remark. Experimentation suggests that this family A n might in fact yield a tight upper bound 
for closed 1-vertex triangulations; see Section [J] for details. 



3.2 Bounded triangulations with many normal surfaces 

The number of vertex normal surfaces in a bounded triangulation has a direct impact on algorithms 
such as unknot recognition and incompressibility testing [19]. Here we give the first explicit expo- 
nential lower bound on the worst-case growth rate of this quantity. The proof is based on a general 
construction principle (Lemma [2]) which, for an arbitrary bounded triangulation Qq satisfying cer- 
tain weak constraints, uses a recursive squaring argument to obtain a family of triangulations {Qk} 
with (/?"■) vertex normal surfaces, where the exponential base (3 is derived from Qq. By choosing 
a suitable starting triangulation Qq, we obtain the explicit base /3 ~ 2.3715 (Corollary [3]). 

For both unknot recognition and incompressibility testing, we can improve the underlying al- 
gorithms by only considering vertex normal discs (vertex normal surfaces that are topologically 
trivial). In Theorem U] we show that this restricted quantity is also worst-case exponential: we 
build a family of triangulations with 0(2 n ) vertex normal discs. 

Lemma 2. Suppose Qq is a bounded triangulation with uq tetrahedra, /o is a boundary face of Qq 
such that not all vertices of /o are identified in Qq, and cq is one of the three normal arc types on 
/o • If there are Qo vertex normal surfaces in Qq that meet /o in at least one arc of type cq but in 
no other normal arc types, then Qq can be extended to a family of triangulations {Qk} in which the 
number of vertex normal surfaces grows at a rate ofQ({3 n ), where (3 = aQ^ n ° +1 \ 

In the proof (for a detailed proof see the appendix), we recursively construct Qk by joining two 
copies of Qk-i to an additional tetrahedron A^ along their faces fk-i (see Figure [2]). Suppose there 
are at-i vertex normal surfaces in Qk-i that meet fk-\ in only arcs of type Ck-i- For each pair 
of such surfaces in the two copies of Qk-i, we can combine these surfaces in a way that extends 
through Afc to meet one of the free faces fk of A& in just one chosen normal arc type Ck, and 
this extension yields a vertex normal surface of Qk- There are ct|_ 1 such pairings, and therefore 
«fc > a 1-i sucri vertex normal surfaces of Qk- This recurrence yields the final growth rate of £~2(/3 n ) 

where /3 = ag . We need the assumption that not all vertices of the boundary face /o are 

identified in Qq to show that each triangulation Qk represents a bounded 3-manifold. 
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Figure 2: Attaching two copies of Gk-i to the tetrahedron A. 
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Figure 3: Face pairing graphs of the binomial family Aq, triangulation G, the path family Vq, the 
bounded family £>4, triangulation £, and the closed family C-j. 

Corollary 3. There is a triangulation Q = Go that is the starting point of a family of bounded 
triangulations {Gk}, & > ; with = 12 • 2 k — 1 tetrahedra and o~(Gk) > 2.3715 nfe . 

We prove this using Lemma [2] by choosing a starting triangulation G = Go with no = 11 
tetrahedra and a choice of face fo and arc type cq with ao = 31 643 corresponding vertex normal 
surfaces. This yields a growth rate of £l(/3 n ) with f3 = 31 643 1 / 12 ~ 2.3715. The face pairing graph 
of G is shown in Figure [3] for a detailed construction see the appendix. 

For our final result, we construct the path triangulation V n from n tetrahedra Aj = z(0123), 
i = 1, ...,n, by joining tetrahedra Aj and Aj+i by the map i(012) i— )■ (i + 1)(013). It can be shown 
that each V n is a bounded triangulation whose underlying 3-manifold is the 3-ball (for details see 
the appendix). 

Theorem 4. For each n > 1, V n has 2 n+1 + ( n+1 K"+ 2 ) g 0(2 n ) vertex normal discs. 

We prove this by obtaining explicit recurrences for the number of vertex normal surfaces with 
different normal arcs based on the matching equations. We have essentially two choices for each 
normal arc, and either a triangle or quadrilateral can be added to the previous triangulation giving 
the 0(2 n ) growth rate. Again, see the appendix for details. 

3.3 Lower bounds for the size of normal coordinates 

Here we give exponential lower bounds on the worst-case algebraic complexity ft(T). Our bounds 
follow a Fibonacci growth rate of ri( [(1 + \/5)/2] n ) ~ f2(1.618) n . Understanding the worst-case 
k(T) is important for improving the time and space complexity of normal surface enumeration 
algorithms due to a better handling of the integer arithmetic involved (see Section [T]) . 

To obtain such lower bounds, we first construct a family of bounded triangulations, each con- 
taining a vertex normal disc with coordinates growing exponentially in the number of tetrahedra. 
We then close the bounded family using a constant number of additional tetrahedra so that the ver- 
tex normal surface with exponential coordinates is preserved. In this way, we are able to construct 
two families of triangulations, bounded and closed, with Fibonacci type growth rates for k(T). 
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Figure 4: Left: construction of the 1-tetrahedron triangulation of LST(1,2, 3). Right: the merid- 
ian disc di which is a 6-gon with 1 normal quad, 2 normal triangles and a maximum of 3 edge 
intersections. 



The key objects of the construction are so-called layered solid tori |18j . These are parameterised 
triangulations of the solid torus: the layered solid torus denoted LST(o, b, a + b) has as its boundary 
a triangulation of the torus with exactly three boundary edges, such that the meridian disc (the 
unique disc of the solid torus meeting the boundary in a non-contractible closed curve) intersects 
the boundary edges in a, b and a + b points. 

Layered solid tori are very common tools when constructing triangulations of a given type of 
3-manifold (see [20, 22, 26] for more about constructing 3- manifolds) . The most prominent example 
of a layered solid torus, the one tetrahedron triangulation of LST(1, 2, 3), is shown in Figure HI 

Theorem 5. There is a family B n of bounded 1-vertex triangulations with n tetrahedra, where 
B n contains a vertex normal disc d n with maximum coordinate F n +i, where F^ denotes the k-th 
Fibonacci number. 

The family B n consists of layered solid tori of type LST(F n+ i, F n+ 2, F n+ 3) and each vertex 
normal surface d n is the corresponding meridian disc. See the appendix for details of the proof. 

The key idea for the construction of a closed family C n of n-tetrahedron triangulations containing 
a vertex normal surface with exponentially growing coordinates is to find a small m-tetrahedron 
triangulation £ with the same boundary as the family of layered solid tori B n acting like a type 
of plug. By this we mean that £ contains a normal surface intersecting d£ in the same way as 
d n intersects dB n for all n > 1. This gives rise to a vertex normal surface in C n with maximum 
coordinate greater than or equal to F n +i. Since m is constant, this gives the same asymptotic lower 
bound on k(T) for closed triangulations as for bounded triangulations. 

Remark (The triangulation £). There is a 4-tetrahedron bounded triangulation £ with d£ = dB m 
for all m which contains two vertex normal surfaces s and t that can be combined into a normal 
surface intersecting d£ in the same pattern as d m . The face pairing graph of £ is shown in Figure 
02 and the intersection of s and t with its boundary d£ is shown in Figure EJ A detailed description 
of £ can be found in the appendix. 

Theorem 6. There is a family C n of closed 1-vertex triangulations with n tetrahedra, n > 5, each 
containing a vertex normal surface with maximum coordinate at least F n _3 if n = 2 mod 3 or at 
least 2 F„_3 otherwise. 




Figure 5: Intersection of s (left) and t (right) with d£ 
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Input size n 


Closed 


Closed and 1-vertex 


Bounded 


1 


4 


3 


3 


2 


17 


12 


17 


3 


81 


63 


156 


4 


577 


433 


2 308 


5 


5184 


3 961 


45 046 


6 


57 753 


43 584 


995 920 


7 


722 765 


538 409 


25 225 447 


8 


9 787509 


7148483 


695134018 


9 


139103 032 


99450 500 


19 933 661871 


10 


2 046 869 999 


1430 396 979 




Total 


n < 10: 2196 546 921 


n < 10: 1537 582 427 


n< 8: 721402 915 
n < 9: 20 655 064 786 



Table 1: The number of 3- manifold triangulations in the census 

As outlined above, we construct C n by gluing B n -^ and £ along their boundary tori. If n = 
2 mod 3, the meridian disc d n _4 glued to a combination of s and t yields a vertex normal projective 
plane or, if n = 0, 1 mod 3, twice d n _4 with a combination of s and t yields a vertex normal sphere; 
the maximum coordinates are then as stated. See the appendix for details. 

We note that the vertex normal sphere from above (in the case n = 0, 1 mod 3) is the only 
non- vertex linking normal sphere in C n . Detecting these normal surface types is one of the key 
tasks in important 3-manifold problems such as prime decomposition. Hence, the family C n is an 
example for a case where, in order to prove the existence of such a normal sphere, dealing with 
exponentially large normal coordinates cannot be avoided. This is a hint towards the conjecture 
that these problems are intrinsically hard to solve using normal surface enumeration methods. 

4 Experimental behaviour 

We turn now to an experimental study of the combinatorial and algebraic complexities of vertex 
normal surfaces. Our experimental data consists of all closed 3-manifold triangulations with n < 
10 tetrahedra, and all bounded 3-manifold triangulations with n < 8 tetrahedra (each appearing 
precisely once up to relabelling). As shown in Tabled! this yields almost 3 billion triangulations 
in total (2196 546 921 + 721402 915 = 2 917949 836). We also extract the ~ 1.5 billion 1-vertex 
triangulations from the closed census for additional study. 

Generating exhaustive censuses of all possible inputs requires sophisticated algorithms and 
significant computational resources. The n < 10-tetrahedron census of all closed triangulations first 
appeared in [4] (which also describes some of the underlying algorithms). The n < 8-tetrahedron 
census of all bounded triangulations is new to this paper. Moreover, we have constructed this 
bounded census for n = 9, with over 20 billion triangulations; however, we only use n < 8 for 
our experiments because the subsequent analysis of normal surfaces for n = 9 remains out of our 
computational reach (for n = 8 this analysis already consumed years of CPU time). 

Table [2] summarises our experimental results, and gives worst-case and average-case measure- 
ments (labelled Max and Aug respectively) for the quantities cr(7~) and k(T) in our various settings. 
Each measurement is taken over the relevant census of triangulations from Table [U The Closed 
and Bounded columns refer to all closed or bounded triangulations respectively; in the Closed 1- 
vertex column we restrict our attention to 1-vertex triangulations of closed manifolds, and in the 
Bounded, discs column we only count vertex normal discs (not all vertex normal surfaces). We 
have also measured the algebraic complexity k(T) in the 1-vertex and discs-only settings, but we 
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Input 


Combinatorial complexity cr(T) 


Algebraic complexity 




size n 


Closed 


Closed 1-vertex 


Bounded 


Bounded, discs 


Closed 


Bounded 




Max 


Avg 


Max 


Avg 


Max 


Avg 


Max 


Avg 


Max 


Avg 


Max 


Avg 


1 


3 


2.0 


2 


1.7 


7 


5.0 


7 


4.0 


1 


1.0 


1 


1.0 


2 


7 


3.9 


4 


3.3 


14 


8.2 


14 


5.2 


2 


1.2 


2 


1.3 


3 


11 


5.5 


8 


4.9 


35 


14.0 


27 


7.0 


3 


1.5 


3 


1.7 


4 


18 


8.8 


16 


7.8 


85 


31.3 


69 


11.6 


4 


1.8 


6 


2.5 


5 


36 


13.3 


32 


12.0 


236 


69.5 


176 


20.2 


7 


2.1 


12 


3.4 


6 


70 


20.8 


64 


18.6 


688 


152.6 


440 


34.8 


10 


2.3 


20 


4.3 


7 


144 


32.2 


128 


28.8 


1943 


376.6 


1109 


61.7 


16 


2.6 


36 


5.8 


8 


291 


50.2 


256 


44.7 


5725 


947.4 


2768 


112.4 


26 


2.9 


65 


7.5 


9 


584 


78.5 


512 


69.4 










42 


3.2 






10 


1175 


123.2 


1024 


108.2 










68 


3.6 






Growth 


2.03" 


1.56" 


2™ 


1.56™ 


2.73" 


2.23™ 


2.45™ 


1.69" 


1.62™ 


1.14™ 


1.81™ 


1.34™ 



Table 2: Experimental worst-case and average-case results 



omit the details due to space constraints; see the appendix for the details. 

The final row of Table [2] gives a "best estimate" of the exponential growth rate of each quantity 
with respect to n (we just list the base of the exponential, ignoring any coefficients or polynomial 
factors). Most growth rates are estimated by linear regressioiJl, though for cases where the worst 
cases matches a known family of triangulations (see below) we give the corresponding known rate. 

We can make some broad observations from Table [2j 

• The average-case scenarios grow at a significantly slower rate than the worst-case scenarios, 
sometimes astonishingly so. This is consistent with past observations in which "typical" tri- 
angulations exhibit significantly smaller complexity properties than expected (see Section [1]). 

• For closed manifolds, 1-vertex triangulations only give a very slight improvement: the worst 
case drops from the 0(17 n / 4 ) ~ 0(2. 03 n ) family described in [3] down to the 0(2 n ) family 
of Theorem [TJ The closeness of these results is surprising, since the theoretical bounds on 
<t(7~) for 1-vertex triangulations are much smaller than the general bounds (see below), and 
algorithms for working with them are often much simpler |18j . 

• Bounded triangulations exhibit higher complexity properties than their closed counterparts. 
For the combinatorial complexity cr(T) this discrepancy is very pronounced — even the average 
case for bounded triangulations is well above the worst case for closed triangulations. This 
is again consistent with past experiences in working with normal surface algorithms |10j . 
Restricting our attention to normal discs (e.g., for unknot recognition) does alleviate this 
problem somewhat. 

Table [3] compares the experimental behaviour of the worst-case cr(T) against its best known 
theoretical lower and upper bounds (here "lower bounds" refers to families of pathological triangu- 
lations with the highest known growth rate of cr(T), such as those constructed in Section [3]). 

Regarding <r(T): the lower bound of 0(17 n / 4 ) ~ 0(2. 03 n ) is known from [3], and the remaining 
three lower bounds are new to this paper (Theorem [H Corollary [3] and Theorem [4} . The first two 
upper bounds are taken from [5]; the final two 0(64 n ) bounds are well known but do not appear 
in the literature (for a proof we refer to the appendix). Regarding k(T): all four lower bounds of 

3 Specifically, we take a weighted linear regression of log a or logK as a function of n. The weights are taken to be 
1, . . . , n, in order to limit the influence of anomalous small cases. 
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Combinatorial complexity a(T) 


Alg 


ebraic complexity 


<T) 




Lower 


Experimental 


Upper 


Lower 


Experimental 


Upper 




bound 


growth 


bound 


bound 


growth 


bound 


Closed 


0(2.03") 




0(14.556") 


0(1.62") 




0(3.17") 


Closed 1-vertex 


0(2") 




0(4.852") 


0(1.62") 




0(3.17") 


Bounded 


0(2.37") 


~ 2.73" 


0(64") 


0(1.62") 


~ 1.82" 


0(31.63") 


Bounded, discs only 


0(2") 


~ 2.45™ 


0(64") 


0(1.62") 




0(31.63") 



Table 3: Summary of worst-case theoretical and experimental results 



0([(1 + \/5)/2] n ) a ©(1.62 n ) are new to this paper (Theorems [5] and [6]), and all four upper bounds 
are taken from [8]. 

Here we see that the experimental growth rates are much closer to the lower bounds than the 
upper bounds; in particular, an arrow (i — ) indicates that the experimental worst-case growth rate 
is identical to the best lower bound (up to a constant factor). This invites the following conjectures: 

Conjecture 7. For closed triangulations, the maximum number of vertex normal surfaces cr(T) 
for any given n grows at an asymptotic rate of 0(17 n / 4 ) ~ B(2.03 n ), and for closed 1-vertex 
triangulations this reduces to 0(2 n ). 

For closed triangulations as well as closed 1-vertex triangulations, the maximum coordinate of 
any vertex normal surface k{T) for any given n grows at an asymptotic rate of B([(l + \/5)/2] n ) ~ 
0(1. 62 n ). For bounded triangulations, the maximum coordinate of any vertex normal disc for any 
given n likewise grows at an asymptotic rate o/0([(l + y/l)/2] n ) ~ 0(1.62"). 

For closed 1-vertex triangulations, our experimental data gives an even stronger result: 

Theorem 8. For closed 1-vertex triangulations, the maximum number of vertex normal surfaces 
cr(7~) for any given n < 10 is precisely 2 n , and is attained by the binomial triangulations A n as 
described in Section \3.1\ 

Conjecture 9. Theorem^ is true for all positive integers n. 

5 Discussion 

As noted in the introduction, although we consider close to 3 billion distinct triangulations, they 
are all relatively small with less than or equal to 10 tetrahedra. Despite this, there are reasons to 
believe that our experimental results might generalise. Because we allow more flexible triangulations 
(not just simplicial complexes), this census contains a rich diversity of 3-manifolds (including 5114 
distinct closed P 2 -irreducible manifolds [U [25]). Moreover, several of the patterns that we see in 
Table [2] (in particular Theorem [8|) are established early on and are remarkably consistent. 

As an exception, in the case of arbitrarily bounded manifolds, a greater number of tetrahedra 
seems to allow a better choice for the starting triangulation Qq in Lemma [2] in order to obtain a 
higher exponential base. 

As seen in Table [31 there is still a long way to go before lower bounds and upper bounds on 
worst-case complexities converge. Not only does this paper produce the first explicit lower bounds in 
several computationally important settings, but it also gives strong experimental evidence that these 
lower bounds are close to (or even exactly) tight. This suggests that it is now the upper bounds that 
require significant improvement, inviting new directions of research with a rich interplay between 
topology, polytopes and complexity theory. 
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Appendix: Additional proofs 



Additional Preliminaries 

In order to fully understand the detailed proofs given in this appendix we need some more termi- 
nology concerning triangulations and normal surface theory. 

From Section [2] we know that normal surfaces within a triangulation can be described by their 
normal coordinates which, due to the matching equations and the non-negativity constraints, lie 
within a high-dimensional polyhedral cone. However, not all of the integer points in this cone give 
rise to a normal surface. Those normal coordinates which in fact do belong to a normal surface are 
called admissible^ Given an admissible vector in this cone, the corresponding normal surface can 
be reconstructed uniquely (up to certain forms of isotopy) . 

The reason why not all normal coordinates are admissible is because there is no way to in- 
sert two quadrilaterals of different types into the same tetrahedron without having a non-empty 
intersection (cf. Figure [I]). Since normal surfaces are defined to be embedded, it follows that for 
every tetrahedron there can be at most one type of normal quadrilateral (though there may be 
many quadrilaterals of this one type). Translated to normal coordinates, this means that for every 
tetrahedron of a triangulation at most one of the three coordinates accounting for the quadrilaterals 
may be non-zero. This condition is called the quadrilateral constraints, and it can be shown that 
the admissible points — that is, the normal coordinates that do correspond to a normal surface — are 
precisely those integer points in the solution cone that satisfy the quadrilateral constraints. 

Inside the cone it is straightforward to define the sum of two normal surfaces as the sum of their 
normal coordinates: if two normal coordinates satisfy the matching equations and non-negativity 
constraints then their sum will also satisfy them. However, not all of these sums represent normal 
surfaces, since they need not satisfy the quadrilateral constraints. Because of this we call two 
normal surfaces compatible if the sum of their normal coordinates is admissible. 

Due to the gluings of a triangulation, many vertices of the tetrahedra may become identified 
together and, as a consequence, a small neighbourhood of a vertex in a triangulation looks like a 
union of cusps of tetrahedra which are themselves smaller tetrahedra. Hence, the outer boundary 
(or frontier) of this neighbourhood is a collection of triangles which all are normal pieces. Therefore 
this outer boundary is a normal surface, which we call a vertex linking normal surface or just a 
vertex link. In a closed triangulation (which represents a closed manifold), all vertex links must be 
spheres (otherwise the neighbourhood of some vertex is not locally M 3 ). In a bounded triangulation 
(which represents a manifold with boundary), all vertex links must be spheres or discs. 

The proof of Theorem Q] 

In this section, we give a full proof of a more detailed version of Theorem [T] from Section 13.11 

Theorem 10. For each n > 1, there are u) vertex normal surfaces in A n of genus k. Therefore 
there are o~(A n ) = 2 n vertex normal surfaces. 

Proof. Consider a single tetrahedron A, with two faces identified as in Section 13.11 Define a, f3 
and 7 to be loops surrounding vertices 0, 1 and 2 respectively. The projective solution space has 
four vertex normal surfaces a, b, c and d. All of these surfaces have a non-empty boundary which 
consists of a combination of a, (3 and 7 and all of them are compatible. Therefore we do not have 
to consider the quadrilateral constraints, which simplifies the rest of the argument. 

4 Some authors define "admissible" to include all non-negative multiples of such points (so admissible vectors need 
not have integer coordinates) In this paper we restrict the notion of admissibility to integer vectors only. 
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Let aj, hi, Cj, dj denote the surfaces and ati, Pi, 7$ the boundary curves in the corresponding 
Aj. Any normal surface in A n can be written as 

n 

^ Ta t &i + T bi bi + T Cl Cj + T d Ai 
i=l 

for constants T ai , 77, . , T Ci , . > for z = 1, ...,n. Following the argument in [3], there are only two 
ways in which a vertex normal surface can meet each Aj and satisfy the matching equations, namely 
Uj = hi + Cj + dj and Vj = aj + dj, both with boundary aj + /3j + 7j. It follows that every normal 
surface in A n can be described by 

n 

^PiUi + gjVj 
i=i 

where each pi, qi > and where e := p± + q± = P2 + 52 = ••• = Pn + q-n is an integer. 

It follows from the above that each normal surface meets the boundary of Aj, i = 1, . . . ,n, in 
e(aj + /3j + 7j). Each normal surface with e > 1 is the sum of two other normal surfaces in A n and 
thus is not a vertex normal surface. On the other hand, if e = 1 it can be shown that the only 
rational linear combinations pi + g$ = 1, i = 1, . . . , n, resulting in a normal surface of A n must be 
Pi = 1 and Qj = or pi = and qi = 1. 

Therefore, the projective solution space has 2™ admissible vertices, corresponding to the 2 n 
combinations of either Pi = 1 and qi = or pi = and qi = 1 at each Aj. The binomial coefficients 
easily follow from a simple Euler characteristic argument since each surface Uj is a disc and each 
Vj is a punctured torus. □ 



Proof of Lemma [2] 

Here we give a full proof of Lemma [2] from Section [3.21 which we restate below. 

Lemma 2. Suppose Qq is a bounded triangulation with uq tetrahedra, /o is a boundary face of Go 
such that not all vertices of /o are identified in Qq, and cq is one of the three normal arc types on 
fo . If there are a® vertex normal surfaces in Qq that meet /o in at least one arc of type cq but in 
no other normal arc types, then Qq can be extended to a family of triangulations {Qk} in which the 
number of vertex normal surfaces grows at a rate ofQ,((3 n ), where (3 = ctQ^ no+1 \ 

Proof. The proof consists of constructing such a family of triangulations. To do so, we take two 
copies of Qq and an additional tetrahedron Ai and join both copies along their faces /o to Ai (see 
Figure E] where k = 1), yielding a new triangulation Q\. Any pair of the olq vertex normal surfaces 
meeting /o i n a c °py of normal arcs of type cq can be combined such that this linear combination 
has a unique extension throughout Ai meeting a third face f\ of Ai in only one of the three normal 
arcs c\. An argument exploiting the uniqueness of this extension shows that any such pairing 
yields a new vertex normal surface and there are a\ = Oq such pairings. The assumption that 
not all vertices of the boundary face /o are identified in Qq is necessary to show that the resulting 
triangulation still represents a bounded manifold (i.e., all vertex links are discs or spheres). 

In a next step, two copies of Q\ are combined together in an analogous way to create Q2, and so 
on. After k steps we have a full binary tree of depth k where the "leaves" of the tree are all copies 
of Qq. In each step the number of tetrahedra is doubled plus one extra "root tetrahedron" is added. 
On the other hand, we at least square the number of vertex normal surfaces. These observations 
put together are sufficient to yield the desired result. 
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Figure 6: Attaching two copies of Qk-i to the root tetrahedron A&. 

More precisely: let Q% be the triangulation after the fc-th step with rik tetrahedra, with a 
designated boundary face fk with a fixed normal arc type Ck of the root tetrahedron A& which will 
be glued to the new root tetrahedron Afc + i, and with at least at vertex normal surfaces meeting 
/it only in copies of cy. (see Figure ED . 

Then we have n^+i = 2n^ + 1 and a^+i > ( a k) 2 - Now let /3 be a real number such that 
> /3 n ° +1 and let us assume that au > /3 rtfc+1 . Since a/c+i > ( a fc) 2 it then follows that 

a k+ i > (P n « +l ) 2 

_ pn k+1 +l 

and thus Gk+i has at least p n k+i+ 1 vertex normal surfaces. To finish the induction argument we 
now simply choose /3 = a^^ no+l \ 

It remains to show that (i) Q}~ is a triangulation of a bounded 3-manifold and (ii) any of the a& 
normal surfaces is in fact a vertex normal surface. 

To show that Qk represents a bounded 3-manifold, we first have to take a closer look at how 
the two copies of Qk-i (denoted by Q\_ x and Q^-i) have to be attached to the two boundary 
triangles of type fk-i (denoted by and f^_i) are joined to A& along an edge e. The gluing 
has to be in a way that the endpoints of e are not previously identified in both fh_i and f%_i, and 
that the normal arcs of type Ck-i of both triangles (denoted by c\_ x and c^_ 1 ) are next to the 
same vertex of A^ which thus has to be an endpoint of e (see Figured]). Note that this is always 
possible if not all vertices of fl_i are identified together and likewise with f"t—\- 

Since Qq = Q is a bounded triangulation, in order to prove that is a bounded triangulation it 
suffices to assume that Qt~i is a bounded manifold and then show that all vertex links of vertices 
of the root tetrahedron A& are still triangulations of the disc after attaching Ql , and Q"i_\ to A&. 

If none of the vertices of ft_i are identified in Gk—u * = 1; 2, none of the vertices of A^ will be 
identified. Thus, all vertex links are either just copies of the vertex links in Qj c _ 1 and Q\_\ with an 
additional triangle attached or two vertex linking discs attached to each other along an additional 
triangle (see Figure [7] on the left). 

If two of the vertices of and f k 2 _ 1 are identified, and if these are the vertices opposite the 
arcs cl_i and c^_ 2 , then one of these vertices in each face must be disjoint from e in A&. In this 
case, three vertices of A& are identified to the same vertex and the vertex link consists of the two 
disjoint vertex links in Ql_i and 5Li joined together along a triangle (the normal piece near the 
endpoint of e opposite of c\_ l and c^_ 2 ), and two additional disjoint triangles which are added at 
opposite ends of the new disc. The result is still a disc (see Figure [7] on the right). The fourth 
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Fi gure 7: Left: no two vertices of fk—\ or 1 identified. Right: two vertices of fk—\ and fk—\ 
are identified. lk(v) denotes the link of vertex v. 

vertex of A/% simply connects two discs as described in the previous situation. It follows that Q/. 
must represent a bounded 3-manifold, and the new face fk likewise does not have all three vertices 
identified. 

If two of the vertices of and are identified, and if these include the vertices inside arcs 
ci , and c\_ 2 , then we must be careful how we attach the two copies of Qk-i to A/%: in one copy 
Gf._i we ensure that the two identified vertices join to the common edge e, and in the other copy 
we only allow one of these vertices to join to e. By a similar argument, this ensures that the 
new vertex links in Qk are discs, and that the new boundary face fk again does not have all three 
vertices identified. 

To prove that all ctk normal surfaces constructed in the fc-th step are in fact vertex normal 
surfaces, consider two normal surfaces ui and U2 meeting fk-i in exactly x\ and X2 normal arcs 
of type Ck-i- Let m be the smallest common multiple of x\ and X2, then ^ ui + ^ U2 can be 
extended to a normal surface in Qk by inserting m matching triangles into A^ (cf. Figure E} . We 
will denote this normal surface by u. Now assume that u can be written as u = As + for two 
normal surfaces s and t of Qk and two rational numbers A and \i. By construction, the restriction 
of this linear combination to Q\_\ or Q\_^ is a multiple of a vertex normal surface. Hence, inside 
or ^fc-i both s and t are rational multiples of u. Moreover, since u inside A& only consists of 
m triangles of a given type, s and t have to be multiples of u inside A^ as well and by the matching 
equations for s and t at fli and fi\ it follows that s and t are multiples of u throughout the 
entire triangulation Qk- Thus, u is a vertex normal surface of Qk- □ 

Construction of triangulation Q 

In Section [3.21 we refer to an 11-tetrahedra triangulation Q as a starting point for the family {Qk} 
of bounded triangulations in Corollary [3l Here we present this 11-tetrahedron triangulation in full. 

The triangulation Q, consisting of n = 11 tetrahedra Aj = i(0123), i = 0, . . . , 10, is given by the 
gluings in Table HI where missing entries denote boundary faces. Face /o is given by /o = 0(012) 
and the normal arc type Co is the one around vertex of fo, hence cq = 0. Using Regina [61 [7] for 
vertex normal surface enumeration yields cr(Q) = 61 526 and ao = 31 643. The face pairing graph 
of Q is shown in Figure [3j 

Proof of Theorem |4] 

Here we give a full proof of Theorem H] from Section 13.2} which we restate below. 
Theorem 4. For each n>\,V n has 2 n+1 + ( TO+1 K n + 2 ) g 0(2") vertex normal discs. 
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A* 


i(012) 


i(013) 


i(023) 


i(123) 







1(012) 


2(021) 




1 


0(013) 


3(012 


4(021) 




2 


0(032) 


5(012) 


6(021) 




3 


1(013) 


7(231) 


7(023) 




4 


1(032) 


8(231 


8(023) 




5 


2(013) 


9(231) 


9(023) 




6 


2(032) 


10(231) 


10(023) 




7 






3(023) 


3(301) 


8 






4(023) 


4(301) 


9 






5(023) 


5(301) 


10 






6(023) 


6(301) 



Table 4: The triangulation Q 




Figure 8: Notation for the normal arcs of the path family. 



Proof. First of all, let V% be the triangulation after the i-th step with i tetrahedra with a designated 
boundary face fi of Aj, and with fixed normal arcs ai, bi and c, on this face which will be glued 
to tetrahedron Aj+i, as shown in Figure Also let on, j3i and 7$ be the number of vertex normal 
surfaces meeting face fi in ai, bi and Cj respectively. Finally, let 5i be the number of vertex normal 
surfaces not meeting face fi. 

First, let us consider the case of a single tetrahedron Ai. There are seven admissible vertices in 
the projective solution space of Ai corresponding to seven vertex normal discs: four single triangles 
around each vertex and three single quadrilaterals. There are two surfaces meeting at each type 
of normal arc in face f\ and one surface without any normal arcs in this face. Therefore a\ = 2, 
P x = 2, 7l = 2 and S 1 = 1. 

Now, let us construct explicit formulae for each of the sequences a» , fy, 7$ and Si . In tetrahedron 
Aj+i = (i + 1)(0123), normal surfaces with normal arcs a.- l+ i cither have a triangle around vertex 
or a quadrilateral separating edge 03 from edge 12, where the notation is as indicated in Figure 
[8l The triangle meets face fi in normal arc aj and the quadrilateral meets fi in normal arc bi. 
Therefore the number of normal surfaces a^+i is exactly a» + fy. A similar argument shows that 
= fa + ai, and since a± = fi\ = 2, it follows that a% = = 2*. 

Normal surfaces with normal arcs q+i on the other hand either have a triangle around vertex 2 
or a quadrilateral separating edge 01 from edge 23 (see Figure [8]). The triangle does not meet fi and 
the quadrilateral meets fi in normal arc Cj. Therefore the number of normal surfaces Ji+i = 7i + 1. 
Therefore since 71 = 2 we have 7^ = i + 1. 

Now consider <5;+i, the number of normal surfaces not meeting face fi+i- These surfaces either 
have a triangle around vertex 3 or do not meet face fi (see Figure [8j) . Since there are exactly <5j 
normal surfaces disjoint to fi and the triangle around vertex 3 meets face fi in a normal arc of type 
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1(012) h> 2(103) 
3 1 . 1 —p. o 

1(023) H- 2(032) 




Figure 9: Building B n in the case n = 2. Left: boundary of d2 in dB<i- Center: the 3rd tetrahedron 
is glued to B2 Right: boundary of d3 in 883 . 



Ci, we have 5i + i =^i + 5%. It follows that 5i = since 5± = 1 and ji = i + 1. 

We note at this point that all vertex normal surfaces in V n have been accounted for, since any 
normal surface that meets the face /j+i in more than one arc can be expressed as a sum of the 
surfaces described above. Therefore the total number of vertex normal surfaces in V n is 

n{n + l) _ on+1 , (n + l)(n + 2) 



a n + Pn + 7n + S n = T + 2 n + (n + 1) + v - ' = 2 n+1 + 



2 2 

It remains to show that each of the above vertex normal surfaces is a disc. At the (i + l)-st step, 
there is only one triangle or one quadrilateral in Aj+i glued to one normal arc in fi. Therefore, 
we either add one vertex, two edges and one face in the case of the triangle or two vertices, three 
edges and one face in the case of the quadrilateral to the surfaces in V%- In both cases the Euler 
characteristic does not change in the (i + l)-st step. Since in tetrahedron Ai all seven surfaces are 
discs, it follows that every surface in V n is a disc. □ 



Proof of Theorem [5] 

In the following, we give a full proof of Theorem [5] from Section 13.31 which we restate below. 

Theorem 5. There is a family B n of bounded 1-vertex triangulations with n tetrahedra, where 
B n contains a vertex normal disc d n with maximum coordinate F n+ i, where F^. denotes the k-th 
Fibonacci number. 

Proof. The family of triangulations B n consists of layered solid tori of type LST(F n+ i, F n+ 2, F n+ 3) 
and each vertex normal surface d n will be the corresponding meridian disc with algebraic complexity 
of at least F n+ i. 

In order to prove this we will start with the base case B\ = LST(1,2,3). Figure [H shows a 
one tetrahedron triangulation of LST(1,2,3) with the meridian disc having maximum coordinate 
F2 = 1. Furthermore, by using Regina (HJ [7] we can check that the meridian disc B n is a vertex 
normal surface of B n . 

Now, let us assume that B n is a n-tetrahedra triangulation of LST(F n+ i, F n+ 2, F ra+ 3) containing 
the meridian disc d n as a vertex normal surface with maximum coordinate F n+ i, and intersecting 
the boundary edges in F n+ i, F n+ 2 and F n+ 3 points (see Figure [9] on the left). 

To construct B n+ \ and d n+ i we now glue the (n + l)-st tetrahedron to LST(F n+ i, F n+ 2> Fn+3) 
such that the boundary edge with F n+ i intersecting points becomes an internal edge (see Figure [9] 
in the center). Note that £> n +i still triangulates a solid torus. There are F n+ 2 parallel normal arcs 
both in the upper left and in the lower right corner (in each case F n+ i arcs coming from triangles 
and F n arcs coming from quadrilaterals of d n ). For each of them we will insert a normal triangle 
into the (n + l)-st tetrahedron of B n +i, intersecting the boundary of B n +i in two arcs and the newly 
inserted boundary edge e in one point. For each of the remaining F n+ i pairs of arcs in the middle 
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of the boundary of B n we insert F n+ i parallel normal quadrilaterals into the (n + l)-st tetrahedron 
of i3 n +i, each intersecting the boundary of £> n +i in two normal arcs and edge e in one point. 

Altogether, we extend the meridian disc of B n to a new normal surface by adding 2 F n+ 2 triangles 
and F n+ i quadrilaterals, hence F n+ 4 normal pieces each intersecting e in one point. The new normal 
surface has two types of normal triangles with each F n+ 2 copies, and thus algebraic complexity of 
at least F n+ 2. Starting from d n , each of these extensions applied one after another turns the 
previous disc into a new disc. As a consequence, the normal surface constructed by attaching all 
F n+ 4 normal pieces to d n is still a disc and will be denoted by d n+ i. The boundary curve of d n+ i 
is still simply closed and non-contractible in the boundary of B n +i, hence, it is the meridian disc 
of B n+ \ intersecting the boundary edges in F„ + 2, F n+ 3 and F n+ 4 points. By construction, B n+ \ 
is a layered solid torus of type LST(F n+ 2, F„ + 3, F„ + 4) with n + 1 tetrahedra (see Figure [9] on the 
right). 

It remains to show that d n+ i is a vertex normal surface. This is equivalent to the statement 
that for any rational combination 

d n +i = As + (1) 

with s and t normal surfaces in B n +i and A, /x 6 Q, the normal surfaces s and t must be rational 
multiples of d n+ \. 

Equation [1] is valid for every subset of normal coordinates in B n +i. In particular, it holds if 
we restrict s, t and d n+ i to the normal coordinates associated to the first n tetrahedra of B n +\. 
However, this restriction applied to d n+ i yields d n which is a vertex normal surface by assumption. 
It follows that s = 6d n and t = ipd n in B n for some rational numbers 6 and if). 

Extending B n to B n+ \, we glue the (n + l)-st tetrahedron onto the two boundary triangles of 
B n , thus adding seven new variables and six new matching constraints to the projective solution 
space of B n in order to enumerate the vertex normal surfaces of B n +%. Following the labelling of 
Figure [9] on the right, we will denote the seven normal coordinates of the (n + l)-st tetrahedron by 

(to, ti, t 2 , h I <?01) <?02, <?03) 

where tx denotes the triangle type isolating vertex X from the rest of the tetrahedron and qxy 
denotes the quadrilateral type isolating edge XY. Taking the normal coordinates from the boundary 
curve of d n into account, this results in the following equations 








= t2 + <?03 







= h + <?03 


Eji- 


h2 


= t + q i 




= h + q i 


En- 


hi 


= t 2 + q 2 




hi 


= h + q 2 



It follows immediately that t 2 = t% = go3 = 0, and qo 2 = F n+ i, qoi = by the quadrilateral 
constraints and thus to = t\ = F n+ 2. Hence, the matching equations assure that d n has a unique 
extension and the same holds for the restricted versions of s and t for some rational multiple of 
the matching equations. It follows that s and t are rational multiples of d n+ i and d n+ i is a vertex 
normal surface. □ 

The triangulation 8 

In Section [3?3] we refer to a 4-tetrahedron triangulation £ acting as a kind of "plug" for the bounded 
family B n in order to obtain the closed family C n with large algebraic complexity. Here we present 
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Figure 10: Left: boundary pattern of s, two disjoint circles. Right: boundary pattern of t 



this 4-tetrahedron triangulation including its vertex normal surfaces s and t in detail. 

The triangulation £ is given by the 4 tetrahedra Aj = i(0123), i = 0, ... ,3, and the following 
table of gluings: 





»(012) 


i(013) 


i(023) 


z(123) 





2(231) 


1(230) 


2(023) 


1(123) 


1 


3(012) 


2(102) 


0(301) 


0(123) 


2 


1(103) 


3(230) 


0(023) 


0(201) 


3 


1(012) 


di 


2(301) 


d 2 



where d\ = 3(013) and d 2 = 3(123) denote the two boundary faces of £ (see Figure [T0|) . The face 
pairing graph of £ is shown in Figure 

Enumeration of all vertex normal surfaces of £ using Regina [6 , 7J yields 13 surfaces, 12 of them 
bounded and one of them closed. However, in the following we will only take a closer look at two of 
the bounded surfaces which will be denoted by s and t. Surface s is a cylinder with two boundary 
components (see Figure [TUl on the left), t is a Mobius strip with one boundary component (see 
Figure [10l on the right). In particular, both surfaces have Euler characteristic 0. 

The surfaces are given by the following standard normal coordinates (cf. proof of Theorem [5] 
for more about the notation of standard normal coordinates): 

s = (0,1,1,0(0,0,0), (0,1,0,0(0,1,0), (0,1,0,0(0,0,1), (0,2,0,2(0,0,0) 
t = (0,0,0,0(1,0,0), (1, 1,0,0 1 0,0,0), (1,1,0,0(0,0,0), (0,1,0,1(0,0,1), 

and boundary patterns 

8s = (0,2,2), (2,0,2) 
dt = (0,2,1), (1,0,2), 

where the first triple counts the normal arcs of boundary face d\ and the second triple the ones of 
boundary face d 2 ■ The triples are ordered starting with the normal arcs around the smallest vertex 
label. 

The surfaces s and t are compatible, meaning that their normal coordinates can be added to 
each other yielding a new normal surface. Since the Euler characteristic is additive under this 
summing operation it follows that all linear combinations of s and t will have Euler characteristic 
0. Moreover, any linear combination of s and t is connected if and only if it is orientable with 
two boundary components or non-orientable with one boundary component. There cannot be any 
interior closed connected components since the only closed vertex normal surface is not compatible 
with t and a scalar multiple of s cannot contain closed connected components. Hence, any connected 
linear combination of s and t is a cylinder or a Mobius strip where the latter is true if and only if 
the linear combination has an odd number of copies of t. 

Keeping this in mind, simple addition of the boundary patterns of s and t shows that 

F n _4 t + - F n _5 s, 
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n yields the boundary pattern of the meridian disc d n _4 from the layered solid torus family 
B n -4 from Theorem [5] above. Moreover, since the boundary patterns of s and t are not multiples of 
each other, this is the only linear combination of s and t with this property. Now, since s contains 
odd coordinates, this boundary pattern bounds a surface if and only if F ra _5 is even, or equivalently 
n = 2 mod 3. In this case F n _4 is odd and since the boundary pattern is connected it bounds a 
Mobius strip which will be denoted by M n . If, on the other hand n = 0, 1 mod 3, the boundary 
pattern 

2 F n _4 t + F n _5 s 

has two connected components and since 2F n _4 is always even, it bounds a cylinder Cyl n . 
Proof of Theorem [6] 

In the following, we give a full proof of Theorem [5] from Section 13.31 which we restate below. 

Theorem 6. There is a family C n of closed 1-vertex triangulations with n tetrahedra, n > 5, each 
containing a vertex normal surface with maximum coordinate at least F n _3 if n = 2 mod 3 or at 
least 2F n _3 otherwise. 

Proof. We will construct the family C n by gluing £>„_4 and £ along their boundary components. 
Then, if n = 2 mod 3, the meridian disc d n _4 glued with M„ yields a vertex normal projective 
plane or, if n = 0, 1 mod 3, twice d n _4 glued with Cyl n yields a vertex normal sphere of C n and the 
maximum normal coordinates are as stated. We will call these surfaces S n . 

It remains to show that S n is a vertex normal surface of C n for all n > 5. To see this, recall that 
d n _4 is a vertex normal surface of £> n _4 for all n > 5. Hence, it suffices to show that d„ (2d„) has 
a unique extension in C n yielding S n if n = 2 mod 3 (n = 0, 1 mod 3) and that S n has the smallest 
integer normal coordinates amongst all rational multiples of S n . 

From the section above we know that there is at least one linear combination of vertex normal 
surfaces in £ realising a valid extension of d„ or 2d„ yielding <S n . Using the classification of 
all vertex normal surfaces of £ we can see that due to conflicting quadrilateral constraints and 
matching equations there are only four vertex normal surfaces in £ which may occur in such a linear 
combination. Due to further compatibility constraints such a linear combination can be shown to 
be a linear combination of s and t and the only closed vertex normal surfaces of £. Following the 
observations made in the section above it follows that there is no other linear combination of vertex 
normal surfaces and hence the extension of d n (2d n ) in £ is unique. 

Altogether it follows that S n is a vertex normal surface. □ 

An upper bound for <x(T) 

The following statement (mentioned in Section^) is well known but does not appear in the literature, 
and so we give the simple proof here. 

Lemma 11. Let T be a triangulation with n tetrahedra. Then o~(T) < 64 n . 

Proof. It is known that each vertex normal surface is uniquely defined by its zero set, i. e., the set 
of normal coordinates which are zero (see [5]). Let (to,ti,t2,ts \ qoi,qo2,Qo3) be the set of normal 
coordinates of a tetrahedron of T ■ Then each of the can be either zero or non-zero which leaves 
us with 2 4 = 16 choices, and at most one of the qoj can be non-zero due to the quadrilateral 
constraints which leaves us with 4 additional choices. All together we have 4 • 2 4 = 64 distinct zero 
sets per tetrahedron which yields the result. □ 
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